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Abstract 



We have done a fast simulation analysis of the H ± signal at LHC in the 1 and 3 prong 
hadronic r-jet channels along with the ti background. The r polarization was effectively 
used to suppress the background in both the channels. Combining this with appropriate 
cuts on the px of the r-jet, the missing E T and the azimuthal angle between them reduces 
the background below the signal level. Consequently one gets a viable H ± signal up to a 
mass range of 600-700 GeV at moderate to large tan /3. 



Introduction 



The Minimal Supersymmetric Standard Model(MSSM) contains a pair of charged Higgs 
bosons H ± along with three neutral ones. While it may be hard to distinguish any of the 
neutral Higgs bosons from that of the standard model, the H^ 1 carries the unambiguous 
hallmark of the MSSM Higgs sector. Therefore, it has an important role in the search of 
MSSM Higgs bosons at the Large Hadron collider (LHC). The leading order QCD process 

gb^t^ + h.c (1) 

gives a sizable production cross section for a heavy at LHC. However its dominant decay 
mode, H^ 1 — > tb, suffers from a large QCD background [1]. A more promising signature 
comes from its leading sub dominant decay mode, 

H± - tV, (2) 

which accounts for a branching fraction of -10% in the moderate to large tan/3( £10) region. 
Moreover one can enhance this signal over the SM background from 

W ± -> r±v r (3) 

by exploiting the opposite polarizations of r i.e. P T = +1 and -1 from the signal(2) and 
background (3) respectively. This was shown to give a viable signature for a heavy H ± boson 
at LHC in its 1-prong hadronic r decay channel [2]. In particular the signal to background 
ratio was shown to be enhanced significantly by requiring the charged prong to carry > 80% 
of the visible r jet energy. 

The work of ref [2] was based on a parton level Monte Carlo simulation for the H ± — > t ± v t 
signal(2) and the W ± — > r^v T background(3), followed by a simple hadronic r decay code 
via 

t -> 7^(12.5%), p ± z/(26%), af 1/(15%). (4) 

The af mode contributes half and half to the 1-prong (n^n !! ) and 3-prong (7r ± 7r ± 7r =F ) 
channels. So the three mesons of eq.(4) account for over 90% of the 1-prong hadronic decay 
branching ratio(BR) of r (~50%) [3]. A more exact analysis was done in [4] following the 
same procedure, where the signal(2) and background(3) were simulated using PYTHIA Monte 
Carlo(MC) event generator [6] along with the fast CMS JET package for detector simulation [7]. 
A similar analysis was also done by members of the ATLAS collaboration [8]. In this paper 
we have investigated the signal(2) and background( 3) in both 1 and 3 prong hadronic decay 
channels of r along the lines of ref. [4] . 

However, we have used the TAUOLA package [9] for hadronic r decay unlike ref. [4], which 
had used the simple decay code of ref. [2] via eq.(4). The two r decay programs give very 
similar results. But the TAUOLA package is more exact, since it includes the small non 
resonant contribution to hadronic r decay. Besides this is the first investigation of this 
process including both 1 and 3 prong hadronic decay channels of r. 
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t Polarization 



It is easy to understand the effect of r polarization(P r ) on its 1-prong hadronic decay via 
the dominant contributions of eq.(4). The center of mass angular distributions of r into it 
or a vector meson v(=p, ai) is simply given in terms of its polarization as 

' ^ -i \ ■ P T COs9) 



T w d cos 9 2 

\l,T 2 ln Ti" L v 



1 dT, - i- 2 ™ 2 



T v d cos 9 m 2 + 2m. 2 . 



;i±P T cos^) (5) 



where L, T denote the longitudinal and transverse polarization states of the vector meson. 
This angle is related to the fraction x of the r lab momentum carried by the meson, i.e the 
(visible) r-jet momentum via 

x = i(i + cos e) + ^(i- COS 9). (6) 

It is clear from eqs.(5) and (6) that the signal (P T = + 1) has a harder r-jet than the 
background (P T — —1) for the it, p L and an contributions; but it is the opposite for p T and 
aw contributions. Now the transverse p and a\ decays favor even sharing of the momentum 
among the decay pions, while the longitudinal p and a\ decays favor uneven distributions, 
where the charged pion carries either very little or most of the momentum. Thus requiring 
the tt^ to carry -80% of the r jet momentum, 

Rl = ll^ > 0.8, (7) 

Vt— jet 

retains about half the longitudinal p along with the pion but very little of the transverse 
contributions. This cut suppresses not only the W — > tv background, but also the fake r 
background to the 1-prong hadronic decay channel from QCD jets. 

The 3-prong hadronic decay accounts for about 15% of r— decay, of which 2/3rd (10%) 
comes from 

t -> 7T ± 7r ± 7r =F I/ (8) 

without any accompanying 7r°. We shall consider only this 3-prong decay channel, which can 
be separated by either matching the tracker momentum of the r-jet with the calorimetric en- 
ergy deposit or by a veto on accompanying n° — > 2j in the electromagnetic calorimeter (EM). 
This effectively suppresses the fake r background from QCD jets while retaining 2/3rd of 
the genuine r events. As mentioned above 3/4th of the 3-prong decay(8) comes from the a\ 
contribution. Thus one can again enhance the an contribution by imposing a cut on the 
fractional r-jet momentum carried by the like-sign pair, i.e 

R 3 = ^ 0.2 - 0.8. (9) 

Vt~ jet 
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We shall see that this cut favors the (P T —+1) signal over the (P T = — 1) background signifi- 
cantly even after including the non-resonant contribution to eq.(8). Note that the like sign 
pion pair in the 3-prong (7r ± 7r ± 7r T ) decay of a\ is analogous to the neutral pion pair in its 
l-prong(7T°7r 7r ± ) decay, i.e. R 3 corresponds to 1-R\ for the a± channel. 

Signal and Background 

We have computed the H ± signal from the leading order QCD process(l) using PYTHIA[6]. 
For simplicity we have used a common renormalization and factorization scale fiR = fiF = s. 
The resulting H ± cross section has been enhanced by a K-factor of 1.5 to account for the 
higher order corrections following ref [10]. Note that the Yukawa coupling of the signal 
process( 1) is estimated using the running quark masses m t (fifi) and m^iiR). Then it is 
followed by the hadronic decay of top, t — > bqq', along with the H — > decay of eq.(2). 
Finally the hadronic decay of r is simulated using TAUOLA [9]. Thus the final state consists 
of at least four hard jets, including the r-jet, along with a large missing- E T (fir)- These 
events are expected to be triggered by a multi-jet trigger along with a higher level r trigger 
with high trigger efficiency of 90% at CMS [4] 1 .The jets and the $t are reconstructed with 
the fast CMS detector response simulation package CMS JET [7]. The program containscthe 
detector resolution and the main cracks and inefficiencies. 

The background was computed from the leading order ti production process using hr = 
\if = s and multiplying the resulting cross-section with the appropriate K-factor 1.3 [11]. 
This is followed by the hadronic decay of one top, t — > bqq', while the other decay via t 
— > bW — > brv. We do not consider the other contributions to the background of eq.(3) coming 
from W+multijet production, since they can be effectively suppressed by the reconstructions 
of the hadronic top mass and b-tagging [4]. Note however that the same polarization cut 
suppresses this W+multijet background as much as the ti background. 

The r identification is based on the narrowness of the r jet. To implement this we define 
a narrow signal cone of size AR$ = 0.1 and an isolation cone of size ARj = 0.4 around the 
calorimetric jet axis, where AR is defined via the azimuthal angle and the pseudo-rapidity 
as 

AR = ^Acj) 2 + ~Arf (10) 

We require 1 or 3 charged tracks inside the signal cone, with \r/\ < 2.5 and p T >3 GeV for the 
hardest track, the former corresponding to the pseudo-rapidity coverage of the tracker. We 
further require that there are no other charged tracks with px >1 GeV inside the isolation 
cone to ensure tracker isolation [4]. This was shown to be adequate to suppress the fake r-jet 
background to the W — > tv events of the CDF experiment in the 1-prong channel but not in 

1 This trigger efficiency falls to 63-64% in full simulation [5]. We have not included it in this work, however, 
since it is based on fast simulation. Its effect will be to rescale the size of signal and background shown in 
the tables as well as figures 1-7 by 2/3rd. Correspondingly, the discovery limits of tan/3 in fig. 8 will increase 
by 25%. 
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the 3-prong channel [12]. In the latter case one has to combine the tracker isolation with the 
requirements of narrow width and small invariant mass(< 1.8 GeV) of the r jet candidate 
using the calorimeter information to suppress the fake r jet background [12, 13, 14], which is 
beyond the scope of the present work. Note however that in the high- pt (~ 100 GeV)range 
of our interest, the 3-prong r-jet is expected to be tagged by the secondary vertex[14]. 
Moreover, we are interested in a sub sample of 3-prong r-jets, without accompanying n°s, 
for which the fake background is relatively small, as we shall see below. 

We ensure the absence of 7r°s by requiring the energy of the 3 charged tracks measured in 
the tracker to match with the calorimetric energy deposit of the r-jet within the calorimetric 
energy resolution, i.e 

AE=\E%*-E%\<10GeV. (11) 

With this cut the r-jet identified by the tracker matches well with the actual r-jet of the 
event generator. Therefore, we shall use the tracker identification of r-jet in both 1 and 3 
prong channels. 

In our simulation, as a basic sets of selection cuts for jet reconstruction, we apply 

jP T > 20 GeV, \rf\ < 4.5 (12) 

for all jets, the latter corresponding to the pseudo-rapidity coverage of the calorimeter. We 
require the minimum separation of 

Ai?>0.5 (13) 

between jets. The missing- Et arising mainly due to the presence of v T accompanied with r 
lepton, is reconstructed using the calorimetric informations, and we set a minimum cut 

E T > 30 GeV. (14) 

The Ifrp is expected to be large for signal as it originates mainly from the massive — > TV 
decay, where as in the case of background it comes from a relatively light W — > tv decay. 

In Fig.l we show the distribution of pt of r-jet (pT~ Jet ) f° r signal in the upper panel 
along with the background from tt in the lower panel. These distributions are subject to 
only the basic selection cuts of eqs.(ll-14) and are normalized for integrated luminosity 
C = 100/6 -1 . The signals are shown for two masses of H ± , m H ± =300 and 600 GeV, for 
tan/3 = 40. Evidently, the higher the mass of H ± , the harder are the r-jets. Notice that 
the signal cross section is several orders of magnitude less than the background even for 
p T T - jet >100 GeV. 

In Fig. 2 we demonstrate the distribution in R 1 ( 7) for 1-prong decay channel of r-jet for 
both signal and background with p^ 3 ^ > 100 GeV. The spillover of the distribution to the 
Ri >1 region reflects the 15-20% uncertainty in the calorimetric energy measurement of the 
r-jet. As discussed above, for the signal process where P T = +1, peaks occur at the two 
ends because of the uneven sharing of energy between pions while in the case of background, 
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for which P T = — 1, a peak occurs at the middle due to the almost equal sharing of energy 
between pions. Note that at the large R\ end of the distributions the dominant contribution 
comes from p L and ir channels. Therefore a cut like R\ >0.8 leads to a very large suppression 
of the background while retaining almost half of the signal events. Note that one would in 
any case require a reasonably hard charged track, corresponding to R\ -0.3, for effective r 
identification and rejection of the QCD jet background [13]. So extending this cut to Ri >0.8 
costs very little to the signal, while it effectively suppresses the P T =—1 background. For 3- 
prong decay of r-jets we present the R 3 distribution in Fig. 3 for both signal and background, 
where R3 is defined by eq.(9). We have shown the distribution for two values of H ± mass as 
before. As mentioned above, the like sign pair of its plays the same role as the 7r° pair in the 
case of 1-prong decay channel of r-jet via a\ i.e., the 7i n n ± channel. Hence, the like sign 
pion pair carries either very little or most of the r-jet energy for the signal(P r = +1), while 
it carries roughly 2/3rd of the r-jet energy for the background(P r = —1). Consequently 
selecting events in the region R% <0.2 and R3 >0.8 helps to suppress the background. We 
found the efficiency due to this selection cut is about 55% for signal and 35% for background. 
Note that there is no constraint on R 3 for r-identification unlike the 1-prong channel. 

Since the H ± production(l) is accompanied by a top quark, it is useful to reconstruct top 
quark mass. We perform the W and top mass reconstructions from hadronic decay modes 
by requiring at least three reconstructed jets in addition to a single r jet. For the W mass 
reconstruction we require the invariant mass of two jets out of them to be 

= mw ± 15 GeV (15) 

and the corresponding W momenta are obtained out of these jet momenta. In case of several 
pairs satisfying this mass band, the pair having invariant mass closest to mw is chosen. The 
top mass reconstruction is performed using this pair and one of the remaining jets and 
demanding 

m r w c j = m t ± 30 GeV (16) 

Notice here that we have not done any kind of b-tagging. However to take care of b-tagging 
we multiply the signal and background cross section by a b-tagging efficiency^) factor 
0.5 [4]. This includes the loss of efficiency due to the reduced coverage of \rj\ <2.5 for this 
jet. Though the tt background has 2 b jets, one still gets efficiency factor 2x0. 5x (l-0.5)=0.5 
by requiring only one b-tag. 

We have also investigated the relative azimuthal opening angle in the transverse plane 
between r jet and the Ifi T vector, which is also connected with the transverse mass via, 

m T = \j2.p T T - jet .$ T (l - cos A0(r - jet,tf T )). (17) 

Since, in the signal process(2) both r-jet and If T are originating from a comparatively 
massive particle, leading to harder r-jet and missing momentum, it is expected that the 
signal distributions in tut and A0(r— jet, I$ T ) will be much broader than the background(3), 
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which can distinguish the signal and background events. In Fig.4 we show the distribution in 
A0(r — jet, tf T ) for the signal for two values of m H =300 GeV and 600 GeV and tan (3 =40 
along with background. These distributions are for 1-prong r-jet events passed by the 
selection cuts: E T T jet >100 GeV, R x >0.8, E T > 100 GeV. In Fig.5 we show the same 
distribution in A0(r— jet, If T ) for r-jets decaying via 3-prong decay modes. In Fig. 6 and Fig.7 
we present the distributions in ttit for 1-prong and 3-prong channels of r-jet respectively. 
Evidently, as expected the background is concentrated at small azimuthal opening angle 
A(f) ~0, while the signal is peaked at the largest opening angle A<p ~ 180°. Likewise the 
background distribution in rnr is restricted to uit < raw, while the signal distribution is 
peaked at much larger tjit and extends all the way up to itlh, which can also be used to 
estimate the H ± mass. Thus either a cut on Acf) or mj> can suppress the background by 
enormous amount without practically any loss of signal events. We will see later that a 
A0 >60° cut brings down the background level to less than the signal. 

In our simulation we have generated 10 6 events each for signal and background processes 
in the r+multijet chanel. In Table 1 we demonstrate the cumulative effect of our selection 
cuts for 1-prong decay channel of r jet, while Table. 2 shows the same for 3-prong decay 
channel of r jet. We present the number of events surviving out of 10 6 generated events 
after each cut. In both the tables we present the results for signal for two sets of values of 
m H ±,tan(3 =(300 GeV,40) and (600 GeV,40) as well as for tt background. First row shows 
the number of hadronically decaying r particles which come from H ± decay of eq. (2) and W 
decay background of eq.(3). It reflects the hadronic r decay branching ratio of ~65%. The 
second row shows the number of events having an identified r-jet, which are subject to the 
basic cuts of eqs. (12-14) and where r-identification is performed using tracker informations 
as discussed above. As expected the corresponding efficiencies for signal are higher than the 
background because of harder r-jet from if ± decay. For the same reason the efficiency after 
the Pr~ Jet > 100 GeV cut, shown in the third row, is about 2.8% for the background where as 
for signal it is 33%(58%) for m H = 300(600) GeV. In the signal events about 3/4 contribution 
comes from the 1-prong decay channel and 1/4 contribution is from 3-prong decay channel 
of r, where as for the background it is 60% from 1-prong channel and 40% from 3-prong 
decay channel of r. The ratio of signal events in the 1 and 3 prong channels agree with the 
respective r branching fractions of 50 and 15%. This shows that the r identification via the 
tracker as described above works quite well for the signal events. On the other hand there is 
a clear excess of 3-prong events in the tt background, showing a large contamination of fake 
r from hadronic jets in this channel, as mentioned earlier. However, they are removed after 
the removal of accompanying 7r°s via the AE cut of eq.(ll). After this cut the 3-prong events 
are about l/5th of the 1-prong events for both the signal and background in agreement with 
the respective r branching fractions of 50 and 10%. 

The next rows show the effects of r polarizations on the signal and background. The 
Ri >0.8 cut retains about 40% of the signal against only 10% of the background. In fact 
the effective loss to the signal is quite small, since the low R\ peak of Fig. 2 would be lost 
anyway due to the requirement of a hard charged track for r identification as mentioned 
earlier. The corresponding cuts of i?3 >0.8 or <0.2 retains about 55% of the signal against 
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35% of the background. It shows the efficacy of r polarization in the extraction of the H ± 
signal in the 1-prong as well as 3-prong r-jet channels. Note also that the measurement of 
the polarization variables R\ and R3 is quite simple, since it only requires measuring the 
momenta of the charged tracks in the tracker. 

The E T >100 GeV cut has an efficiency of 60(90)% for a m H ± =300(600) GeV Higgs 
signal and 40% for the ti background. This is followed by the W and top mass cuts which 
have a combined efficiency 20% for the signal as well as background. Its main utility is in 
suppressing the W+multijet background, as mentioned earlier. Finally, the A0 > 60° cut 
suppresses the background by a factor of 50, with very little loss to the signal. 

The last but one rows show the signal and background cross-sections in the r+multijet 
channel, representing the 10 6 generated events. It corresponds to cr tH ± x Bi(H ± — > tu) x2/3 
for the signal and 2 a t i x 1/9 x 2/3 for the background. The last rows show the signal and 
background cross-sections remaining after all the cuts, which includes a b-tagging efficiency 
factor of 0.5. We see that at this stage the background is reduced to less than the signal size. 
It can be still further suppressed via the transverse mass distribution without any loss of 
signal. The tut distribution can also be used to estimate m#, as mentioned earlier. Thus the 
discovery limit is primarily controlled by the size of the signal shown in the last row. Fig. 8 
shows the discovery limit corresponding to 25 signal events for 1-prong and (l+3)-prong 
channel of r-jet. This is shown separately for C = and 100/ft" 1 expected from the 

low and high luminosity runs of LHC separately. It shows a promising discovery potential for 
over the mass range upto 600-700 GeV at moderate to large tan /3. However it should be 
noted here that a full simulation for the high luminosity run has not been completed yet [14]. 
It would require higher trigger threshold of r-jets, which would move up the discovery limits 
at low ttih to slightly higher values of tan (3. 

Summary 

We have investigated the signal at LHC in the 1 and 3 prong hadronic r-jet channels 
along with the tt background. The signal and background processes were generated using 
the PYTHIA event generator interfaced with TAU0LA for the 1 and 3 prong hadronic r decays. 
The detector response was simulated using the fast simulation packages CMSJET. We have 
shown that the opposite polarization of r from the signal (H^ — > tv) and background 
(W ± — > tv) processes can be effectively used to suppress the background with respect to 
the signal in both 1 and 3-prong r-jet channels. The signal was also found to have much 
harder distributions than the background in the azimuthal opening angle(A0) between the 
r-jet and the missing Et as well as in the transverse mass (%) of the two. Combining these 
distinctive features with those of r polarization we could effectively suppress the background 
to below the signal size. Thus the H ± discovery potential of LHC in the 1 and 3 prong r-jet 
channels is primarily determined by the signal cross section in these channels. We find a 
promising H ± signal at LHC over the mass range of several hundred GeV at moderate to 
large tan (3. 
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Table 1: Number of events after each set of cuts for two sets of (M H ,tan (3) values for 
the signal and the tt background process. Number of r+multjet events generated in each 
case is 10 6 . Last two rows show the production cross sections in the r+multijet channel and 
cross section after multiplying with efficiency factors including the b-tagging efficicency for 
the 1-prong hadronic r decay channel. 



Cuts 


ttih, tan/5 


m,H, tan (5 


Bg 




(300,40) 


(600,40) 


tt 


No. of had.r decay event 


640531 


641346 


641009 


Identified tau jets 


448288 


465556 


229622 


E T T 3et >100 GeV 


212491 


370562 


17875 


1 prong decay 


158865 


281428 


10182 


R lw >0.8 


57240 


103470 


1046 


E T > 100 GeV 


32441 


89587 


410 


Number of jets >3 


15083 


41783 


308 


W mass rec, rrijj = raw ± 15GeV 


9861 


27091 


147 


Top mass rec,m,jw — m>t ± 30GeV 


6376 


17339 


87 


A0(r - jets, tf T ) > 60 u 


5154 


16394 


2 


axBR (pb) 


0.431 


0.045 


73 


Cross section x efficiency xe&(fb) 


1.1 


0.37 


.15 



Table 2: Same as in Table 1, but for r jets in 3-prong decay channel. 



Cuts 


mu, tan (5 


mu, tan (3 


Bg 




(300,40) 


(600,40) 


tt 


No. of had.r decay event 


640531 


641346 


641009 


Identified tau jets 


448288 


465556 


229622 


E T T ~ jet >100 GeV 


212491 


370562 


17875 


3 prong decay 


53626 


89134 


7683 


AE < 10 GeV 


32886 


54901 


2610 


R 3n <0.4 or > 0.8 


18159 


29714 


858 


E T > 100 GeV 


10456 


26173 


269 


Number of jets >3 


4854 


12161 


206 


W mass rec, rrijj = mw ± 15GeV 


3138 


7886 


110 


Top mass rec, rrijw — n^t =t 30GeV 


2010 


5073 


60 


A0(r - jets, E T ) > 60° 


1676 


4881 


1 


crxBr(pb) 


.431 


0.045 


73 


Cross section x efficiency xe b (fb) 


0.36 


0.11 


0.07 
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Figure 1: The number of events are shown against px of r-jets for signal and background 
processes for integrated luminosity C =100/& 1 and are subject to p^ 3 ^ >20 GeV and 
\rj T -j et \ <2.5 cuts. The mass of charged Higgs and tan/5 are set to m H ± =300 and 600 GeV 
and tan (3 =40. 
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Figure 2: The number of events are plotted against the fraction of r-jet momentum carried 
by the charged prong R\ for the 1-prong decay channel of r-jet for signal and background 
processes. Both the distributions are subject to PT~ jet >100 GeV and \rj T -j et \ <2.5 cuts. The 
masses of charged Higgs and tan (3 are same as in Fig.l. 
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Figure 3: Same as in Fig. 2, but against the fraction of r-jet momentum carried by the like 
sign pair(i? 3 ) for the 3-prong r-jet channel. 
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Figure 4: The number of events are shown against the opening azimuthal angle A0(r — 
jets, If T ) for signal and background for 1-prong decay channel of r-jets. These are subject 
to p T jT jet >100 GeV, >0.8 and tf T > 100 GeV cuts. 
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Figure 5: Same as in Fig.4, but for 3-prong decay channel of r-jets. 
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Figure 6: The number of events are shown against transeverse mass tut for signal and 
background for 1-prong decay channel of r-jets. These are subject to p^~ jet >100 GeV, 
R\ >0.8 and I$ T > 100 GeV. The masses of charged Higgs and tan (3 are same as in previous 
figures. 
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Figure 7: Same as in Fig.5, but for 3-prong decay channel of r jets. 
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Figure 8: Discovery limits of charged Higgs are shown as functions of tan/3 for integrated 
luminosity C = 30 fb -1 (solid lines) and C = 100 fb -1 (dashed lines) for H ± — > rv where 
r decays hadronically in both 1-prong and 3-prong channels. In each case the contribution 
from 1-prong channel only are shown by upper lines where as contribution from the combined 
1 and 3-prong channels are shown by lower lines. 
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